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Abstract 

The quantum hydro-dynamic model for charged particle systems is ex- 
tended to the cases of non zero magnetic fields. In this way, quantum 
corrections to magnetohydrodynamics are obtained starting from the 
quantum hydrodynamical model with magnetic fields. The quantum 
magnetohydrodynamics model is analyzed in the infinite conductivity 
limit. The conditions for equilibrium in ideal quantum magnetohydro- 
dynamics are established. Translationally invariant exact equilibrium 
solutions are obtained in the case of the ideal quantum magnetohy- 
drodynamic model. 

P ACS numbers: 52.30.Cv, 52.55.-s, 05.60.Gg 



1 Introduction 

There has been an accrued interest on quantum plasmas, motivated by appli- 
cations in ultra small electronic devices pQ, dense astrophysical plasmas 
and laser plasmas [SJ. Recent developments involves quantum corrections to 
Bernstein- Greene-Kruskal equilibria quantum beam instabilities [?]-[TU]. 
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quantum ion-acoustic waves jllj . quantum corrections to the Zakharov equa- 
tions modifications on Debye screening for quantum plasmas with 
magnetic fields |14j, quantum drift waves [T3], quantum surface waves [TE], 
quantum plasma echoes [T7|, the expansion of a quantum electron gas into 
vacuum [TH] and the quantum Landau damping [TH]. In addition, quan- 
tum methods have been used for the treatment of classical plasma problems 

[Ulini- 
One possible approach to charged particle systems where quantum effects 
are relevant is furnished by quantum hydrodynamics models. In fact, hydro- 
dynamic formulations have appeared in the early days of quantum mechanics 
|22| . More recently, the quantum hydrodynamics model for semiconductors 
has been introduced to handle questions like negative differential resistance 
as well as resonant tunneling phenomena in micro-electronic devices [23J- 
|25j . The derivation and application of the quantum hydrodynamics model 
for charged particle systems is the subject of a series of recent works [26J- 
In classical plasmas physics, fluid models are ubiquitous, with their 
applications ranging from astrophysics to controlled nuclear fusion j^l IS] ■ 
In particular, magnetohydrodynamics provides one of the most useful fluid 
models, focusing on the global properties of the plasma. The purpose of this 
work is to obtain a quantum counterpart of magnetohydrodynamics, start- 
ing from the quantum hydrodynamics model for charged particle systems. 
This provides another place to study the way quantum physics can modify 
classical plasma physics. However, it should be noted that the quantum hy- 
drodynamic model for charged particle systems was build for non magnetized 
systems only. To obtain a quantum modified magnetohydrodynamics, this 
work also offer the appropriated extension of the quantum hydrodynamics 
model to the cases of non zero magnetic field. 

The paper is organized as follows. In Section 2, the equations of quantum 
hydrodynamics are obtained, now allowing for the presence of magnetic fields. 
The approach for this is based on a Wigner equation with non zero vector 
potentials. Defining macroscopic quantities like charge density and current 
through moments of the Wigner function, we arrive at the desired quantum 
fluid model. In Section 3, we repeat the well known steps for the derivation 
of magnetohydrodynamics, now including the quantum corrections present 
in the quantum hydrodynamic model. This produces a quantum magneto- 
hydrodynamics set of equations. In Section 4, a simplified set of quantum 
magnetohydrodynamics is derived, yielding a quantum version of the general- 
ized Ohm's law. In addition, the infinity conductivity case is shown to imply 
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an ideal quantum magnetohydrodynamic model. In this ideal case, there is 
the presence of quantum corrections modifying the transport of momentum 
and the equation for the electric field. Section 5 studies the influence of the 
quantum terms on the equilibrium solutions. Exact solutions are found for 
translational invariance. Section 6 is devoted to the conclusions. 



2 Quantum Hydrodynamics in the Presence 
of Magnetic Fields 

For completeness, we begin with the derivation of the Wigner-Maxwell sys- 
tem providing a kinetic description for quantum plasmas in the presence of 
electromagnetic fields. For notational simplicity, we first consider a quantum 
hydrodynamics model for non zero magnetic fields in the case of a single 
species plasma. Extension to multi-species plasmas is then straightforward. 
Our starting point is a statistical mixture with N states described by the 
wave functions ip a = ip a (r, t), each with probability p a , with a — 1 ... N. Of 
course, p a > and J2a=iPa — 1- The wave functions obey the Schrodinger 
equation, 

1 fhb 
— {-ihV - qA) 2 i, a + q<j>i/> a = ih^ . (1) 

Here we consider charge carriers of mass m and charge q, subjected to pos- 
sibly self-consistent scalar and vector potentials <fi = <f>(r,t) and A = A(r, t) 
respectively. For convenience in some calculations, we assume the Coulomb 
gauge, V • A = 0. 

From the statistical mixture, we construct the Wigner function / = 
/(r, p, t) defined as usual from 

1 N r s s 

f(*^t) = ——J2p a ^;(r + -)eTi( r --). (2) 

V / a=l 

After some long but simple calculations involving the Schrodinger equations 
for each ip a and the choice of the Coulomb gauge, we arrive at the following 
integro-differential equation for the Wigner function, 

f + £ 'V/= (3) 

at m 

'* 1 fdsdp'e li2 ^[^T+^- ( p(v-^]f(r,p',t) + 



h(2nhf J J * L ^ v 2' ^ v 2' 
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V • / / efe dp' [A(r + J) - A(r - J)] /(r, p', t) 



2m(27rh) 3 J J 2 2 

" P • / / * e • [A(r + -) - A(r - -)] /(r, p', t) . 

All macroscopic quantities like charge and current densities can be found 
taking appropriated moments of the Wigner function. This is analogous 
to classical kinetic theory, where charge and current densities are obtained 
from moments of the one-particle distribution function. Alternatively, we 
could have started from the complete many body wave function, defined 
a many body Wigner function and then obtained a quantum Bogoliubov- 
Born-Green-Kirkwood-Yvon hierarchy. With some closure hypothesis, in this 
way we arrive at a integro-differential equation for the one-particle Wigner 
function, which has to be supplemented by Maxwell equations. This is the 
Wigner-Maxwell system, which plays, in quantum physics, the same role the 
Vlasov-Maxwell system plays in classical physics. When the vector potential 
is zero, it reproduces the well known Wigner-Poisson system |45[ l4*p r j . In 
addition, in the formal classical limit when h — ► 0, the Wigner equation Q 
goes to the Vlasov equation, 

^ + v.V/4(E + vxB)io, (4) 
ot m av 

where v = (p — qA)/m, E = — V0 — dA/dt and B = V x A. However, 
notice that a initially positive definite Wigner function can evolve in such a 
way it becomes negative in some regions of phase space. Hence, it can not 
be considered as a true probability function. Nevertheless, all macroscopic 
quantities like charge, mass and current densities can be obtained from the 
Wigner function through appropriated moments. 

Equation (jSJ) coupled to Maxwell equations provides a self-consistent ki- 
netic description for a quantum plasma. As long as we know, it has been first 
obtained, with a different notation, in the work |17j. It has been rediscovered 
in [3H] , in the case of homogeneous magnetic fields. Wigner functions appro- 
priated to non zero magnetic fields have also been discussed, for instance, 
in jM]-[ST], without the derivation of an evolution equation for the Wigner 
function alone. More recently, a different transport equation for Wigner func- 
tions appropriated to non zeros magnetic field and spin has been obtained in 
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[52] . The starting point of this latter development, however, is the Pauli and 
not the Schrodinger equation as here and j^Jj. Finally, relativistic models 
for self-consistent charged particle systems with spin can be found in [53 a . 

Most of the works dealing with quantum charged particle systems prefer 
to work with the wave functions and not directly with the Wigner function, 
as in The impressive form of (j3J) seems to support this approach. In- 
deed, probably (J3j) can be directly useful only in the linear or homogeneous 
magnetic field cases. This justifies the introduction of alternative descrip- 
tions. At the coast of the loss of some information about kinetic phenemena 
like Landau damping, we can simplify our model adopting a formal hydro- 
dynamic formulation. Define the fluid density 

n = Jdpf, (5) 

the fluid velocity 

u = — [dp(p-qA)f (6) 

mn J 

and the pressure dyad 

P = At / dp (p - qA) <g> (p - qA) f - mi <g) u . (7) 

m z J 

We could proceed to higher order moments of the Wigner function, but ©- 
(J7J) are sufficient if we do not want to offer a detailed description of energy 
transport. 

Taking the appropriated moments of the Wigner equation Q and using 
the definitions ©-(0), we arrive at the following quantum hydrodynamic 
model, 

V • (mi) = , (8) 

u • Vu = ——V • P + — (E + u x B) . (9) 

n m 

Equations (jEJ)-© does not show in an obvious way any quantum effects, 
since h is not explicitly present there. To found the hidden quantum effects, 
we follow mainly the style of references > E±] and [33] , but now allowing 
for magnetic fields. In the definition (J2J) of the Wigner function, consider the 
decomposition 

ip* = V^e iSa/h , (10) 



+ 



on 

~dt 
du 

at + 
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for real n a = n a (r, t) and S a = S a (r, t). Evaluating, the integral for the pres- 
sure dyad, we get a decomposition in terms of "classical" P c and "quantum" 
P^ contributions, 

P = P C + P Q , (11) 



where 



N 



m Pa n a(u a - u) <g> (u a - U) + (12) 
a=l 

N 

+ m$> a n a «-u°)®«-u°), 

0=1 

h 2 n 

P Q = V®Vlnn. (13) 

4m 

In the definitions of classical pressure dyad P c , we considered the kinetic 
fluid velocity associated to the wave function ip a , 

U Q = , (14) 

m 

and the kinetic fluid velocity associated to the statistical mixture, 

JV 

u a . (15) 

a=l ^ 

In a similar way, the second term at the right hand side of equation (|12|) is 
constructed in terms of u°, the osmotic fluid velocity associated to the wave 
function ip a , 

U *=2^—> (16) 
and u°, the osmotic fluid velocity associated to the statistical mixture, 

N 



a=l 



u o = yP^a u o (17) 

n 



We also observe that in terms of the fluid density n a of the state a the density 
n of the statistical mixture is given by 

N 

n = Yl p* n » ■ ( 18 ) 

a=l 
6 



Notice that P a faithful classical pressure dyad, since it comes from dis- 
persion of the velocities, vanishing for a pure state. Indeed, the classical 
pressure dyad is the sum of a kinetic part, arising from the dispersion of 
the kinetic velocities, and a osmotic part, arising from the dispersion of the 
osmotic velocities. However, P c is not strictly classical, since it contains h 
through the osmotic velocities. In a sense, however, it is "classical" , since it 
comes from statistical dispersion of the velocities. 

In most cases, it suffices to take some equation of state for P c . For sim- 
plicity, from now on we assume a diagonal, isotropic form = SijP, where 
P = P(n) is a suitable equation of state. Certainly, strong magnetic fields 
have to be treated more carefully, since they are associated to anisotropic 
pressure dyads. However, since we are mainly interested on the role of the 
quantum effects, we disregard such possibility here. 

Now inserting the preceding results for the pressure dyad into the mo- 
mentum transport equation Q, we obtain the suggestive equation 



The equation of continuity (jSJ) and the force equation (fTTJj) constitute our 
quantum hydrodynamic model for magnetized systems. All the quantum 
effects are contained in the last term of the equation the so called Bohm 
potential. In comparison with standard fluid models for charged particle 
systems, the Bohm potential is the only quantum contribution, and the rest 
of the paper is devoted to study its consequences for magnetohydrodynamics. 



3 Quantum Magnetohydrodynamics Model 



The equations from the last Section were written for a single species charged 
particle system. Now we generalize to a two species system. Consider elec- 
trons with fluid density n e , fluid velocity u e , charge — e, mass m e and pressure 
P e . In an analogous fashion, consider ions with fluid density rii, fluid velocity 
Uj, charge e, mass m ; and pressure p. Proceeding as before, now starting 
from the Wigner equations for electrons and ions, we get the following bipolar 
quantum fluid model, 




du 



ran m 



— VP + — (E + ux B) + 





(20) 
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^ + V-(n iUi ) = 0, (21) 
du P _ VP P e 



+ u e Vu e = (E + u e xB) + 



dt ' vn e n e m. 

n 2 _ /v\/rn 



— -^ + Ui-Vu, = - + — (E + Ui xB) + 

2 /V72 



+ ^vf^-^j-i/^-Ue). (23) 

In the equations (|22M23|) . we have added some often used phenomenological 
terms to take into account for the momentum transport by collisions. The 
coefficients v ei and v ie are called collision frequencies for momentum transfer 
between electrons and ions |%3*1 HI] . For quasineutral plasmas, global mo- 
mentum conservation in collisions imply m e v ei = mji/ ie , so that v ie <C v ei 
when the ions are much more massive than electrons O • 

Equations (j2T?|) - (f2^|) have to be supplemented by Maxwell equations, 

VE = (24) 
VB = 0, (25) 

f)~D 

VxE = -_, (26) 

dE 

VxB = fi J + ^ e Q —, (27) 

where the charge and current densities are given respectively by 

p = e (rij - n e ) , J = e (71*11* - n e u e ) . (28) 



Equations (|2()H28j) constitute our complete quantum hydrodynamic model, 
allowing for magnetic fields. When B = 0, it goes to the well known quantum 
hydrodynamic model for bipolar charged particle systems. 

Several possibilities of study are open starting from (|2UH28p . Here we are 
interested in obtaining equations analogous to the classical magnetohydrody- 
namic equations. In some places, for the sake of clarity and to point exactly 
for the new contributions of quantum nature, we repeat some well known 
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+ ^n^j + ^ v ht"j' (32) 



n = P I + mW ^ (Ue - u,) ® (u e - ui) , (33) 

Pm 



steps in the derivation of classical magnetohydrodynamics. To proceed in 
this direction, define the global mass density 

p m = m e n e + miUi (29) 

and the global fluid velocity 

m e n e u e + miUiUi 
U = — . (30) 

m e n e + miiii 

With these definitions and proceeding like in any plasma physics book jJHJ 
, we obtain the following equations for p m and U, 

^ + V-(p m U) = 0, (31) 

Prni-QjT + U-VU) = -V-n + pE + JxB + 



for 



where P = P e + Pi and where I is the identity matrix. In equations 
the electronic and ionic densities are defined in terms of the mass and charge 
densities according to 

n e = ; [Pm P) , Ui = ■ [Pm H P) ■ (34) 

rrii + m e e rrii + m e e 

We can simplify (J32|) considerably assuming, as usual, quasi-neutrality (p = 
so that n e = P e = Pi = P/2 and neglecting m e in comparison to rrii 
whenever possible. In addition, disregarding the last term at the right hand 
side of (jSni), we obtain 

^ + D .VU=-ivP + IjxB + J^v(^V (35) 
dt p m p m 2m e m, \ v /p^ J 

Under the same assumptions and following the standard derivation of mag- 
netohydrodynamics ^3 E] , we obtain the following equation for the current 
J, 

m e rriidJ mS/P „ „ m; „ h 2 „ /V 2 ^/p^\ 1 , . 

E+UxB -JxB V — v£ll\ j ? (36) 



p m e 2 dt p m e p m e 2em e \ ^fp~^ J a 
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where a = p m e 2 / l {m e m i i' ei ) is the longitudinal electrical conductivity. Equa- 
tion (|36|) is the quantum version of the generalized Ohm's law (SJ El] . The 
continuity equation (j3"Tj) . the force equation (J3l)j) . the quantum version of the 
generalized Ohm's law (jHBjl . an equation of state for P, plus Maxwell equa- 
tions, provides a full system of quantum magnetohydrodynamic equations. 
However, it is probably still complicated and in the next section we propose 
some approximations in the same spirit of those of classical magnetohydro- 
dynamics. 



4 Simplified and Ideal Quantum Magnetohy- 
drodynamic Equations 

Usually O Hlji the left-hand side of the equation (JSTIj) is neglected in the 
cases of slowly varying processes and small pressures. Also, for slowly varying 
and high conductivity problems , the displacement current can be neglected 
in Ampere's law. Finally, we assume an equation of state appropriated for 
adiabatic processes. This provides a complete system of simplified quantum 
magnetohydrodynamic equations, which we collect here for convenience, 



dp m 

dt 



+ V-(p m U)=0, (37) 



^ + U ■ VU = — —VP + —J x B H — — V(^—^=^) , (38) 

at p m Pm 2m e mi ^fp^ 

VP = K 2 Vp m , (39) 

f)T2 

V x E = -— , (40) 

V x B = /i J, (41) 

J = a[E + UxB-^JxB-— V(^)]. (42) 

p m e 2em e ^/p~^ 

In equation ()39j). V s is the adiabatic speed of sound of the fluid. Gauss law 
can be regarded as the initial condition for Faraday's law. Also notice that 
the Hall term J x B at (}4*2"j) is often neglected in magnetohydrodynamics. 

Inserting (J3"9~j) into (|3*SJ). we are left with a system of 13 equations for 13 
unknowns, namely, p m and the components of U, J, B and E. This is our 
quantum magnetohydrodynamics model. In comparison to classical magne- 
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tohydrodynamics, the difference of the present model rests on the presence 
of two quantum corrections, the last terms at equations (|55j) and (|32J) . 

In the ideal magnetohydrodynamics approximation, we assume an infinite 
conductivity and neglect the Hall force at (}4*2|) . This provides the following 
ideal quantum magnetohydrodynamics model, 

E = -UxB + (43) 
2em e ^/p^ 

Pm {— + U-VU) = -VP + — (V x B) x B + 

Ot fJ,Q 

2m e rrii Jp~^ 



V x (U x B) , (45) 



dB 

~dt 

supplemented by the continuity equation (|3*7)l and the equation of state 

Taking into account (fHHJ) . equations (|44H45j) plus the continuity equation 
provides a system of 7 equations for 7 unknowns, namely, p m and the compo- 
nents of U and B. This is our ideal quantum magnetohydrodynamics model. 
In comparison to classical ideal magnetohydrodynamics, the difference of the 
present model rests on the presence of a quantum correction, the last term 
at equation (jUJ). Interestingly, taking the curl of (J4*3~|) makes disappear one 
of the quantum correction terms present in the non ideal quantum magne- 
tohydrodynamics. This leads to a dynamo equation (J45|) identical to that 
of classical magnetohydrodynamics. Consequently, for infinite conductivity 
the magnetic field lines are still frozen to the fluid, even allowing for the 
quantum corrections proposed here. In fact, even for finite conductivity, the 
diffusion of magnetic field lines is described by the same diffusion equation 
as that of classical magnetohydrodynamics. This comes from the fact that 
the quantum correction disappear after taking the curl of both sides of ()42|) . 
neglecting the Hall term and assuming a constant a as usual. However, a 
further quantum correction on the electric field still survives through (|4*3*j) . 

In order to obtain a deeper understanding of the importance of quantum 
effects, we propose the following rescaling for our ideal quantum magnetohy- 
drodynamic equations, 

p m = pjpo, IJ = U/V A , B = B/B , 
f = SIp/Va, i=Qit, (46) 
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where po and B are the equilibrium mass density and magnetic field. In 
addition, Va = (-Bq/(//oPo)) 1,/2 is the Alfven velocity and fij = eB /rrii is the 
ion cyclotron velocity. We justify the chosen rescaling in the following way. 
In magnetohydrodynamics, the Alfen velocity provides a natural velocity 
scale. Also, since we deal with low frequency problems, f^ -1 is a reasonable 
candidate for a natural time scale. These velocity and time scales induces 
the length scale Va/^i, as shown in (|JBJ). 

Applying the rescaling ()46j) to the ideal quantum magnetohydrodynamic 
equations, we obtain the following non dimensional model, 



dp r , 



+ V-(p m U)=0, (47) 



dt 

d\J V 2 - 

p m (^ + U-VU) = -^Vp m + (VxB)xB + 

+ ^V(^), (48) 



OB 

~dt 



V x (U x B) , (49) 



where 



H = Jgl T/2 (50) 

is a non dimensional parameter measuring the relevance of quantum effects. 
Numerically, using M.K.S. units, we have H = 3.42 x 10~ 30 uq/Bq, where no 
is the ambient particle density. While for ordinary plasmas H is negligible, 
for dense astrophysical plasmas with n about 10 29 — 10 34 m~ 3 , H 

can be of order unity or more. Hence, in dense astrophysical plasmas like the 
atmosphere of neutron stars or the interior of massive white dwarfs, quantum 
corrections to magnetohydrodynamics can be of experimental importance. 
Similar comments apply to our non ideal quantum magnetohydrodynamics 
model. However, even for moderate H quantum effects can be negligible if 
the density is slowly varying in comparison with some typical length scale, 
due to the presence of a third order derivative at the Bohm potential. This 
is in the same spirit of the Thomas-Fermi approximation. 
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5 Quantum Ideal Magnet ost at ic Equilibrium 



There is a myriad of developments based on classical magnetohydrodynam- 
ics (linear and nonlinear waves, dynamo theory and so on) and we shall not 
attempt to reproduce all the quantum counterparts of these subjects in the 
framework of our model. We will be restricted to just one subject, namely 
the construction of exact equilibria for ideal quantum magnet ohydro dynam- 
ics, with no attempt to study the important question of the stability of the 
equilibria. 

Assuming that U = and that all quantities are time- independent, the 
ideal quantum magnetohydrodynamic equations (|43H45j) becomes 



h 2 V 2 fo~ 

E = -^V(-^), (51) 



2em e yfp ri 

VP = i(VxB)xB + ^V(^). (52) 
Ho 2m e m; J~p^ 



According to (|51J1. the equilibrium solutions of ideal quantum magnetohydro- 
dynamics are not electric field free any longer. In addition, equation ()52)1 has 
an quantum correction that invalidate the classical magnetic surface equation 
for B • VB = 0, namely P + B 2 /(2fx ) = cte. 

Equation (J52J) together with an equation of state is the key for the search 
of equilibrium solutions. We will try to follow, as long as possible, the strat- 
egy of reference for classical magnetostatic equilibria. Inspired by well 
known classical solutions jSH], assume a translationally invariant solution of 
the form 

P = P(r,cp), p m = p m (r,<p), (53) 
B = VA(r,ip) x z + B z (r,ip)z, (54) 

using cylindrical coordinates and where A = A(r,ip) and B z = B z (r,(p) as 
well as the pressure and the mass density are functions of (r, if) only. 

Substituting the proposal (j53H54j) into we get, for the radial and 

azimuthal components of this equation, 

V (P + JS.) = _! VA V 2 A + V(^l=) , (55) 



2/i Po 2m e m. 
while, for the z component, the result is 

d(B z ,A) 

— r = 0. (56) 

d(r, <p) 
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In ()56|) and in what follows, we used the definition of Jacobian, 



8B Z dB z 

dr dip 

d A d A 

dr dip 



(57) 



d(B z ,A) 
d(r, <p) 

From (JHEJ), we obtain 

B Z = B Z {A). (58) 

Taking into account (|55j) and the fact that B z is a function of A, it follows 
that 

d(P,A) = h 2 Pm d(V 2 ^ ,A) 
d(r, f ) 2m e rrii d(r, <p) 

In the classical limit ft — > 0, the right hand of (j^Sj) vanishes, implying just 
the functional relationship P = P(A). In the present work, we still postulate 



so that, from (|59j). we have 



P = P(A) , (60) 



V2 ^=F(A), (61) 



where F = F(A) is an arbitrary function. 

The last equation is a distinctive feature of ideal quantum magnetohy- 
drodynamic equilibrium. Indeed, ()61|) would not be necessary if h = in 
()59|) . Hence, even if ft is not present in (|6ip. this equation has a quantum 
nature, with important implications in what follows. The reason why ft does 
not appear in (pTjl is that it factor at the right hand side of (|59|) . 

From (jnnj) and some subjacent equation of state, P = P(p m ), we deduce 

Vp^ = G(A) , (62) 
for some function G = G(A). Plugging this into (jBTjl . the result is 

|vM + |(Vi) 2 = F(A), (63) 

where the prime denotes derivation with respect to A. 
Coming back to (jB3j) . we obtain 

V 2 A = ^oh^(A) + G 2 F'(A)] , (64) 
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where we have defined 

K = K{A) = P{A) + . (65) 

2/io 

Recapitulating, we have three four functions of A to be stipulated, namely 
F, G, K and P. However, A satisfy two different equations, (jo^|) and ()64|) . 
Once A is found, all other quantities (pressure, mass density, electromagnetic 
field) comes as consequences. 

A reasonable choice is to take G as a linear function of A, since then (I63J) 
becomes linear in the derivatives. Hence, let 

G = k x A + k 2 , fci ^ , (66) 

for numerical constants k\ and A^- We take k\ ^ since fci = would imply 
F = 0, making disappear the quantum correction at ()64|) . With the choice 
(jnHJ), the couple (|63H64j) becomes 

V 2 A = }-(k x A + k 2 )F(A), (67) 

V 2 A = fiQ [-K'(A) + -^—(k 1 A + k 2 ) 2 F'(A)]. (68) 

2m e mi 

The right hand sides of (|6T|) and (ffiT^) should coincide, implying 

K\A) = (k,A + k 2 ) 2 F\A) - -\- (k x A + k 2 ) F(A) . (69) 

2m e rrii /io^i 

The last equation define K up to an unimportant numerical constant. 

Equation (JBTj) is the key equation for our translationally invariant mag- 
netostatic equilibria. For a given F(A) and solving (J67|) for A, all other 
quantities follows for a known equation of state. Indeed, knowing A we can 
construct the radial and azimuthal components of the magnetic field through 
(I54J) and the mass density from ()62j). From the mass density and the equa- 
tion of state, we obtain the pressure P. Proceeding, equation (jp^j) yields 
K(A) and then the z component of the magnetic field through (J65j) . Finally, 
the electric field follows from (|51|) and the current density from the curl of 
the magnetic field. The free ingredients to be chosen to construct explicitly 
the exact solution are the function F(A) and the equation of state, and the 
numerical constants k\ and k 2 . Other possibilities can be explored if we do 
not restrict to linear G(A) functions as in (J6TI|) . but then A will not satisfy 
an linear in the derivatives equation. 
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5.1 An Explicit Exact Solution 

An interesting case of explicit solution is provided by the choice 

F(A) = klB " { } ~f k) e~ 2kA ^ , (70) 

K\A + k 2 

where B is an arbitrary constant magnetic field, k is an arbitrary constant 
with dimensions of an inverse length and < e < 1. With the choice ([70)1 . 
the equation (}6T|) traduces into the Liouville equation, 

V 2 A= (l-e 2 )B ke- 2kA/Bo , (71) 

which admits the exact cat eye solution 

A = — — ln[cosh(fcr cosip) +e cos(fcr sin <p)] . (72) 
k 

All other relevant quantities follows from this exact solution following the 
recipe just stated. The mass density, from (Jo^j) . is 

k\B 

Pm = [ — ; — ln(cosh(fcr cosy?) + e cos(kr sin tp)) + k 2 ] 2 , (73) 
while the radial and azimuthal components of the magnetic field follows from 

(El, 

g B [sin ip sinh(A;r cos <p) + e cos ip sin(/cr sin <p)] 

[cosh(kr cos if) + e cos(kr simp)] 
g B [cos f sinh(/cr cos^) — e sin <p sm(kr sirup)] 

19 [cosh(/cr cos <p) + e cos(Ar sin ip)] 

Assuming an adiabatic equation of state, P = V s p m , we get, from (fHHjl . 

B 2 Z = B 2 - 2/i V 2 (hA + k 2 ) 2 + (76) 
+ (1 - e 2 ) k 2 e~ 2A [1 + poh (h + k 2 )h 2 /m + /i h 2 k\ A/m] , (77) 

with A given by the cat eye solution (fT2~j) . If desired, the electric field and the 
current density can then be calculated via ()51|) and Ampere's law respectively. 
In figure 1, we show the contour plot of the function A given by (|72|). while 
in figure 2 we show the corresponding mass density. The parameters chosen 
were B — 1, k — 1, e — 0.9, k\ = 1 and k 2 = 0. These graphics shows 
coherent, periodic patterns resembling quantum periodic solutions arising in 
other quantum plasma systems Similar graphics can be easily obtained 
for the electromagnetic field and other macroscopic quantities derivable from 
the cat eye solution (J72J). 
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6 Conclusion 



In this work, we have obtained a quantum version of magnetohydrodynam- 
ics starting from a quantum hydrodynamics model with nonzero magnetic 
fields. In view of its simplicity, this magnetic quantum hydrodynamics model 
seems to be an attractive alternative to the Wigner magnetic equation of Sec- 
tion 2. The infinite conductivity approximation leads to an ideal quantum 
magnetohydrodynamics. For very dense plasmas and not to strong magnetic 
fields, the quantum corrections to magnetohydrodynamics can be relevant, 
as apparent from the parameter H derived in Section 4. Under a number of 
suitable assumptions, we have derived some exact translationally invariant 
quantum ideal magnetostatic solutions. More general quantum ideal mag- 
netostatic equilibria can be conjectured, in particular for axially symmetric 
situations. In addition, we have left a full investigation of linear waves to 
future works. 
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